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ABSTRACT 
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^ I Context. The majority of magnetic chemically peculiar (mCP) stars exhibit periodic light, radio, spectroscopic and spectropolarimetric 

C/2 , variations that can be adequately explained by the model of a rigidly rotating main-sequence star with persistent surface structures. 

CU Vir and V901 Ori belong among these few mCP stars whose rotation periods vary on timescales of decades. 
Aims. We aim to study the stability of the periods in CU Vir and V901 Ori using all accessible observational data containing phase 
information. 

Methods. We collected all available relevant archived observations supplemented with our new measurements of these stars and 
analysed the period variations of the stars using a novel method that allows for the combination of data of diverse sorts. 
Results. We found that the shapes of their phase curves were constant during the last several decades, while the periods were changing. 
At the same time, both stars exhibit alternating intervals of rotational braking and acceleration. The rotation period of CU Vir was 
gradually shortening until the year 1968, when it reached its local minimum of 0.52067198 d. The period then started increasing, 
reaching its local maximum of 0.5207163 d in the year 2005. Since that time the rotation has begun to accelerate again. We also found 
much smaller period changes in CU Vir on a timescale of several years. The rotation period of V901 Ori was increasing for the past 
quarter-century, reaching a maximum of 1.538771 d in the year 2003, when the rotation period began to decrease. 
Conclusions. We propose that dynamical interactions between a thin, outer magnetically-confined envelope, braked by the stellar 
wind, with an inner, faster rotating stellar body is able to explain the observed rotational variability. A theoretically unexpected 
alternating variability of rotation periods in these stars would remove the spin-down time paradox and brings a new insight into 
structure and evolution of magnetic upper-main-sequence stars. 



Key words, stars: chemically peculiar - stars: variables - stars: individual: CU Vir, V901 Ori - stars: rotation 

1 . Introduction centrated into large, persistent spot regions. As the star rotates, 

periodic variations in the brightness, spectrum, and magnetic 

In addition to mass, all stars inherit a fraction of the angular ^^j^ observed. Combining both present and archival obser- 

momentum of the mother cloud out of which they condensed. ^^^^^^^ ^(^p ^^^^^ collected over the past several decades, we 

Consequently, all stars rotate. Stars spend the prevailing part of reconstruct their rotational evolution with unprecedented ac- 

their lives as main-sequence (MS) objects fusing hydrogen in curacy 
their cores. During the MS phase, the angular momentum of stars 

without mass-loss is conserved. This implies that stellar rotation Careful period analyses of photometric observations of 
period changes should occur on the timescale of the moment dozens of the best-monitored mCP stars listed in the On-line 
of inertia variations, which is of the order of 10^ - lO"* years catalogue of photometric observations of magnetic chemically 
(Meynet & Maeder 2000). How do we test this? peculiar stars (Mikulasek et al. 2007) have been peiformed in 
The magnetic chemically peculiar (mCP) stars, which have Mikulasek et al. (201 1). They confirm our theoretical expecta- 
an abnormal surface chemical composition, are the most suit- tions that the rotation periods and light curves of most upper MS 
able test beds for studying the rotational evolution in upper MS stars are constant on timescales of decades. A few mCP stars dis- 
stars. The overabundant elements in their atmospheres are con- play secular cyclic changes in the shape of their light curves (see 
e.g. Ziziiovsky et al. 2000), which can be attributed to the pre- 
sent/ offprint requests to: Zdenek Mikulasek, cession of magnetically distorted stars (Pyper & Adelman 2004). 
e-mail: mikulas@physics. muni . CZ However, there is also a small subgroup of hot mCP stars that 
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have stable light curve shapes and spectroscopic variability, but 
exhibit variable rotation periods (Pyper et al. 1998; Mikulasek et 
al. 2008; Townsend et al. 2010). Constancy of their light curves 
on the scale of decades excludes precession as the cause of the 
observed period changes (for details see Mikulasek et al. 2008). 

This letter concentrates on the period analysis of the best- 
monitored mCP stars - CU Virginis and V901 Orionis. 

2. The stars 

2.1. CU Virginis 

CUVir = HD 124224 = HR 5313 is one of the most enig- 
matic stars in the upper MS. This very fast rotating silicon mCP 
star has a moderate, nearly dipolar magnetic field (Borra & 
Landstreet 1980; Pyper et al. 1998) with the pole strength of 
Bp = 3.0 kG tilted towards the rotational axis by j8 = 74°, the 
axis inclination being i ^ 43° (Trigilio et al. 2000). CUVir is the 
only known MS star that shows variable radio emission, resem- 
bling the radio lighthouse of pulsars (Trigilio et al. 2008, 2011; 
Ravi et al. 2010), it also displays variations in light and spectral 
lines of He i, Si ii, H i, and other ions. The nature of its variabil- 
ity is minutely studied in Krticka et al. (2011). CUVir belongs 
to the most frequently studied mCP stars. 

Occasional rapid increases in its rotation period have been 
reported and discussed. Pyper et al. (1997, 1998) discovered an 
abrupt increase of the period from 0?5206778 to O?52070854 that 
occurred approximately in 1984 and Pyper & Adelman (2004) 
discussed two possible scenarios of the explanation of the ob- 
served O-C diagram, namely a continually changing period or 
two constant periods. After the year 1998 Trigilio et al. (2008, 
201 1) observed another increase in the period of radio pulses of 
AP = 1.12s with respect to the period determined by Pyper et al. 
(1998). However, the 2010 measurements indicate some period 
decrease. The rate of the deceleration can be evaluated using the 
spin-down time, r, defined as r = P/P, where Pit) is the instan- 
taneous rotation period at the time t and P is the mean rate of 
the rotational deceleration. The paradox of CU Vir is according 
to Mikulasek et al. (2011) that its spin down time, r ~ 6 x 10^ 
years, is more than two orders of magnitude shorter than the esti- 
mated age of the star - 9 x 10^ yr (Kochukhov & Bagnulo 2006). 

2.2. V901 Orionis 

V901 Ori = HD 37776 is a young hot star (B2IV) residing in 
the emission nebula IC 432, with a global, extraordinarily strong 
(Bs « 20 kG), and unusually complex magnetic field (Thompson 
& Landstreet 1985; Kochukhov et al. 2011). The observed mod- 
erate light variations are caused by the spots of overabundant 
silicon and helium (Krticka et al. 2007). 

Three decades of precise photometric and spectroscopic 
monitoring have revealed a continuous rotational deceleration 
(Mikulasek et al. 2008), increasing the period of about l'?5387 
by a remarkable 18 s! Ruling out (a) the light-time effect in a bi- 
nary star, (b) the precession of the star's rotational axis, and (c) 
evolutionary efi'ects as possible causes of the period change, we 
interpreted the deceleration in terms of the rotational braking of 
the outer stellar layers as caused by the angular momentum loss 
in the stellar magnetosphere. 

However, this cannot explain the discrepancy between the 
spin-down time, T = 2.5x10^ yr, and the star's age of one milUon 
years or older (Mikulasek et al. 2008, 2010, 2011; Kochukhov 
& Bagnulo 2006). The interpretation of the rotation period evo- 
lution by simple angular motion loss was also questioned by 
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Fig. 1. Long-term variations of the rotation period and O-C values for 
the magnetic chemically peculiar star CU Vir. Left: The model curve of 
long-term rotational period changes AP of CUVir in seconds with re- 
spect to the mean period = O!52069415. The formal accuracy of the 
polynomial fit is comparable with the thickness of the fitted line. Right: 
The difference of the observed and calculated times of zero phase ac- 
cording to a linear ephemeris in days. Each point represents the aver- 
age of 498 consecutive individual measurements; the weights of these 
means are indicated by their areas. 

the negative value of the second derivative of the period: P = 
-29(13) x 10"'^ d"', which indicated that braking could change 
into acceleration (Mikulasek et al. 2008). Nevertheless, the sig- 
nificance of this conclusion was fairly low. 

3. Observations 

For the period analyses of CU Vir we obtained some new pho- 
tometric or spectroscopic data, specifically: 1) JK derived 210 
individual spectrophotometric magnitudes in 21 bands (201 - 
299 nm) by processing 10 lUE spectra, in 1979 January - March 
(for details see Krticka et al. 201 1); 2) GWH obtained 374 pre- 
cise BV measurements at Fairborn Observatory Arizona, USA, 
in 2010 February - 2011 June; 3) JL obtained 38 V measure- 
ments using his own microtelescope at the SAAO, South Africa 
and Brno, CR, in 2010 April - June ; 4) JJ obtained 25 1 vb ob- 
servations at Suhora Observatory, Poland, in 2011 March; and 
JJ-i-JL obtained 536 f/BV measurements using the 0.5m reflector 
of the SAAO, South Africa, in 201 1 May. 

We collected a total of 8965 individual measurements of 
CUVir obtained between 1949 and 2011 (62 years or 43662 
revolutions of the star) including 8270 photometric and spec- 
trophotometric measurements in photometric bands from 200 to 
753 imi as weU as spectroscopic, spectropolarimetric and radio- 
metric observations. The diverse observational data originated 
from 39 sources (see the list in Appendix B). 

The V901 Ori data hsted in Mikulasek et al. (2008) were en- 
hanced with 98 precise BV measurements obtained by GWH at 
Fairborn Observatory Arizona, USA from September 2009 to 
February 20 1 1 . We used a total of 26 1 1 individual measurements 
obtained from 1976 to 201 1 from 14 sources. In addition to 2409 
photometric measurements, the data include 202 measurements 
of equivalent widths of eleven selected He: lines. 

4. The models and their results 

The method of data processing, outlined in the Appendix A, as- 
sumes that phase curves of all monitored quantities are constant. 
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while the 'instant period' of changes P(t) is variable in time t. It 
is advantageous to introduce a 'phase function' &(t, b), which is 
a continuous monotonic function of time t. The fractional part of 
it corresponds to the common phase (p, the integer part is the so- 
called epoch (E) and b is a set of g\, free parameters describing a 
model of the period development. Quantities P(t) and §{t, b) are 
bound by the simple equaUty P(t) -!/■&, where the dot refers to 
the first time derivative of the quantity. Because we did not know 
the true nature of observed period changes of P(t) we chose stan- 
dard low-orders polynomials for their description. 

4.1. CU Virginis 

The period changes P{t) (see Fig. 1) can be adequately approxi- 
mated by a cubic parabola with the origin at the time / = Tq, 



P{t) = Po+a{30-40^); = it-To)/n, 



(1) 



where A and U are parameters expressing the amplitude and 
timescale of the observed period changes, and 6* is a time-like 
function. The instant period P{t) reaches its local extrema Pq+A 
at 0i_2 = +5, where Pq is the period at the time of 1(0 = 0) = Tq. 
The value 77 then equals the duration of the rotation deceleration 
epoch that took place in the time interval f 1 2 = 7"o + ^n. 

Applying the equahty = I /Pit), we obtain the following 
approximation for the phase function in the form of the fourth- 
order polynomial of time 



Po 



t-Mp 
Po ' 



B : 



^, (2) 
7^0 



where ??o is the phase function for a linear ephemeris with the 
origin at Mo = 2 446 730.4447 and the basic period Po- 

Using the above formulated phenomenological model, we 
derive Po = O?52069415(8), 77 = 13260(70) d = 36.29(19) 
yr. To = 2446636(24) or 1986.56(7), A = 1.915(10) s, and 
B = 0.5643(29) d. The period P{t) reached its local minimum in 
the year 1968.4, Pmi„ = 0^52067198(7), and its local maximum 
in 2004.7: Pmax = O?52071631(18). The rotational deceleration 
rate reached its maximum P = 0.158 syr~^ in the year 1986.6. 
The zero phase times of CU Vir can be evaluated using the rela- 
tion JD{k) s Mo + Pok + B i^Ok' - 0k^), where k is an integer 
and 0k ^{Mo+Pok- To)/n. 

The basic model fits the observed long-term period changes 
very well and enables us to predict the zero phase times with an 
accuracy of 0.001 d. However, the analysis of the residuals from 
the accepted polynomial model reveals an additional variation. 
Between 1988 and 1998, for which there is excellent photomet- 
ric coverage (Adelman et al. 1992; Pyper et al. 1998), we find a 
small modulation of the period on a timescale of several years 
(Fig. 2). 

4.2. V901 Ononis 

For the phase function of V901 Ori we adopted polynomials 
in the form of the third-order Taylor-expansion, also used in 
Mikulasek et al. (2008): 

&(t, b) = ^0 - jh^il - iPo P&l; &0 = (f - Mo)/Po; 

P(t) = 1 - Po ( 1 + Po&o + jPoP&l) ' (3) 

where &o is the phase function for the Unear ephemeris with the 
origin at Mo. Po and Po are the instantaneous period and its first 
derivative at Mo. We assumed that the second derivative of the 
period, P, is constant throughout the interval of observations. 




Fig. 2. Medium-term variations in the rotational period of CUVir. The 
modulation of the long-term rotational changes in the observed time 
of the light minimum minus the calculated time of light minimum is 
expressed in days. Owing to uneven sampling within the observing 
seasons, we were not able to detect rotational modulation on shorter 
timescales. 




1990 2000 2010 
Years 



1980 1990 2000 
Years 



Fig. 3. Long-term variations in the rotational period and O-C diagram 
of V901 Ori. Left: Variation in the rotation period AP in seconds with 
respect to the mean period P(, = 1*1538756(3), derived from 2611 in- 
dividual measurements. The 1-cr deflections from the quadratic fit are 
denoted by dashed lines. Points with error bars correspond to virtual 
period deflections for appropriately selected groups of consecutive indi- 
vidual data items. Right: The difference of the observed zero phase time 
and the time calculated according to the linear ephemeris in fractions of 
a day. Each point represents an average of 435 consecutive measure- 
ments; the weights of these means are indicated by their areas. 



We put the origin. Mo = 2449 967.969(5), at a bright- 
ness maximum (in 1999). Our fitted parameter values include 
Pq = l'?538756(3), the instantaneous period at time Mq, and 
Po = 11.3(0.7) X 10-'' = 0.356(23) syr-i and P = -4.4(6) X 
lO^^^d"' - -1.38(19) X lO^^syr"^, the first and second time 
derivatives. The negative value of P is now established with 
7-cr certainty. Accordingly, the star reached its longest period, 
p„^^ = l'?5387709(15), in 2002.8 ± 1.1, and is now accelerat- 
ing again (see Fig. 3). The zero-phase timing can be calculated 
by the relation JD(k) ^ Mo + Po{k+ \Pq + ^PoPk^ where 
k is an integer. There is no indication of any short-term period 
modulation analogous to that seen in CU Vir. 



5. Discussion 

We were able to reveal unexpected alternating lengthenings and 
shortenings in the periods of CUVir and V901 Ori owing to 
simultaneous processing of all available observational material 
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obtained during many decades, including our own observations 
made in 2007-1 1. The period evolutions of both stars were mod- 
elled by polynomials of low degrees. In the case of CU Vir we 
also tested a fit with a harmonic function (Mikulasek et al. 201 1) 
and a Gaussian function corresponding to the transient nature of 
the phenomenon. The only differences among them are notice- 
able at the very begirming of our observations, which shghtly 
favour the latter two alternatives. 

The possible harmonic-like changes of the observed peri- 
ods of both stars together with the strict constancy of the phase 
curves can in principle be explained by the light time effect 
caused by an invisible companion orbiting the star in a nearly 
circular orbit. However, this orbital motion would be revealed 
by considerable radial velocity variations ARV. In the case of 
CUVir the observed change in the period AP of 3.8 s corre- 
sponds to a variation ARV = cAP/P = 25 kms"\ while the 
observed change in V901 Ori period of AP ^ 20 s predicts 
a ARV variation of even 45 kms M However, no such long- 
term RV variations synchronised with period changes have been 
found (Pyper et al. 1997; Mikulasek et al. 2008). Additionally, 
the light time effect cannot explain the complex medium-term 
period variations observed in CU Vir (see Sec. 4. 1, Fig. 2). 

We conclude that the observed period variations are caused 
by uneven rotation of surface layers of the stars (Stgpieh 1998; 
Mikulasek et al. 2008, 2011). However, the physical explana- 
tion of the observed rotation period variations is not straightfor- 
ward. Rotational braking by angular momentum (AM) loss via 
a magnetised stellar wind (Mikulasek et al. 2008) cannot be the 
sole cause of the period variations because we detected intervals 
when the rotational period decreases. Consequently, the AM loss 
may affect only the outer stellar envelope that is hardened by the 
global magnetic field, leaving the faster rotating core unaffected. 
Intervals of rotational deceleration may then alternate with inter- 
vals of angular momentum exchange between the slowly rotating 
envelope and the faster rotating core. The latter manifests itself 
by the rotational acceleration of the outer envelope. 

However, this explanation is theoretically challenging. Given 
the intermediate age of CU Vir, the detected rotation variations 
are likely not connected with the transient behaviour that is 
connected with settling into the MS stage. Therefore, we can 
expect a stable inner configuration in both stars. As follows 
from Flowers & Ruderman (1977) and Braithwaite & Nordlund 
(2006), the magnetic fields confined to the outer stellar layers 
only are unstable. On the other hand, a magnetic field that pen- 
etrates deep into a star causes nearly uniform rotation (Moss et 
al. 1990; Charbonneau & MacGregor 1992; Maeder & Meynet 
2003). 

The possible role of stellar winds in the mechanism of the 
rotational instability is supported by the fact that stronger period 
variations are detected in V901 Ori, which is more luminous than 
CU Vir and, consequently, has a stronger stellar wind (Krticka & 
Kubat 2001). Moreover, this could help to explain the fact that 
the younger star V901 Ori rotates more slowly than CU Vir. 

All these considerations should be the subject of additional 
physical modelling. Nevertheless, the aim of this letter is more 
unpretentious - we only want to draw the attention to an inter- 
esting property of the rotation of some upper MS stars. 



6. Conclusions 

Our study investigated the nature of the rotational period vari- 
ations in two well-observed mCP stars CUVir and V901Ori. 
Contrary to the results of previous studies, we show that the 



period changes are not monotonic - intervals of rotational de- 
celeration alternate with intervals of rotational acceleration, all 
on the timescale of several decades. These results explain the 
spin-down time paradox of these stars, according to which the 
spin-down time was significantly shorter than the age of these 
stars. On the other hand, this unexpected behaviour of two fairly 
dissimilar mCP stars poses a strong challenge for any theoretical 
models. 
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Appendix A: The outline of the method 

The techniques used to analyse the data are based on the 
rigourous application of non-Unear, weighted, least-squares 
methods used simultaneously for all relevant data that contain 
phase information. Our technique does not use an O-C diagram 
as an intermediate stage of data processing; O-C diagrams are 
used only as a visual check on the adequacy of the models. 
Let us assume that all observed phase curves of a star are ade- 
quately described by the unique general model function F{§, a), 
described here by g,, parameters contained in a parameter vector 
a, a = (ai, ...,aj, ...,agj . In our computation we assume that 
the form of all the phase variations is constant and that the time 
variability of the observed quantities are given by a phase func- 
tion §(t, b), which is a monotonic function of time t. The frac- 
tional part of it corresponds to the common phase, the integer 
part is the so-called epoch (E). We can express the phase func- 
tion by means of a simple model quantified by gb parameters b, 
b - {b\, ...,bk, —,bg^). The instantaneous period is given simply 
by the equality P{t) - 

For a realistic modelling of phase variations for all data types 
for a star, we need ga free parameters for the description of the 
model function F{§, a) and gh free parameters for the description 
of the phase function -t^it, b) . The computation of the free param- 
eters was iterative under the basic condition that the weighted 
sum 5(a,b) of the quadrates of the difference Ay, of the ob- 
served value y, and its model prediction is minimal (w, being 
the individual weight of the i-\h measurement). 

n 

Ayi=yi-F{-»i)\ S =Yj^y}wi; 55=0; ^ (A.l) 
i=i 

1=1 I— 1 

We obtain here g - ga+ gb equations of g unknown param- 
eters. The weights of individual measurements w,- are inversely 
proportional to their expected uncertainty. The system is non- 
linear; we have to determine the parameters iteratively. With a 
good initial estimate of the parameter vectors a and b, the itera- 
tions converge very quickly. Usually we need only several tens 
of iterations to complete the iteration procedure. 

A.1. Virtual O-C diagrams. Evoiution of periods. 

The short-term modulation of the phase function was analysed 
by means of the residuals of the observed data Ay,, creating in- 
dividual values of the phase shifts expressed in days (0-C)j with 
adapted individual weight Wj for each observed datum and aver- 
ages of the phase shifts defined for arbitrarily selected groups of 
measurements (O-C), or deflection of the mean period from the 
instant model period APi,(tk): 

(O-C), = -Pa,) Ay, (^j ; W,' = (^) w,; (A.2) 

_ (O-C), Wj _ (O-C), Wj 

(O-C),- ^„^^^^ ; KP,-—^^-^——. 

Computations of (O-C), and AP^ followed after the model 
parameters were found; consequently they had no influence on 
the model solution. They were used only to visualise of the so- 
lution. Similarly, we can compute virtual 'observed' values of 



the instant period from a group of observations to generate the 
model curves in our figures. 

Relation A.l can also be used to determine reliably zero- 
phase times for selected groups of observations O,. These quan- 
tities depend only marginally on the chosen model of the phase 
function §{t, b). Therefore we can use these O, in the process of 
the phase function modelUng. 

Appendix B: Brief specification of CU Vir data 



Table B.l. Data with phase information about CU Vir used in the period 
analysis. 



season/s 


Data type 


N 


Reference 


1949-52 


EWHei 


34 


Deutsch (1952) 


1955 


UBV 


162 


Hardie (1958) 


1964-6 


UBV 


702 


Abuladze (1968) 


1964 


EWSin 


40 


Peterson (1996) 


1964 


EWHei 


9 


Peterson (1996) 


1965-6 


RVHy 


16 


Abt&Snowden(1973) 


1966 


B 


33 


Hardie et al. (1967) 


1967,76 


EWHei 


37 


Hardorp & Megessier (1977) 


1967,76 


EW Sin 


19 


Hardorp & Megessier (1977) 


1968-9 


EWSin 


124 


Krivosheina et al. (1980) 


1968-9 


EWHy 


21 


Krivosheina et al. (1980) 


1968 


UBV 


287 


Blanco & Catalano (1971) 


1972-81 


spf mag 


179 


Pyper & Adelman (1985) 


1972 


UBV 


29 


Winzer(1974) 


1974 


P 


95 


Weiss etal. (1976) 


1974 


uvbyW 


460 


Weiss et al. (1976) 


1975 


spf AO II 


54 


Molnar&Wu(1978) 


1976-8 




14 


Borra & Landstreet (1980) 


1977 


EWHei 


69 


Pedersen (1978) 


1979 


spf lUE 


210 


Krtickaetal. (2011) 


1980-3 


uvby 


89 


Pyper & Adelman (1985) 


1980-3 




23 


Pyper & Adelman (1985) 


1983 


EWHei 


9 


Hiesberger et al. (1995) 


1987-9 


UBV 


1012 


Adelman etal. (1992) 


1988 


P 


25 


Musielok et al. (1990) 


1990-3 


Bj VjHp 


288 


ESA(1998) 


1991-7 


uvby 


3635 


Pyper etal. (1998) 


1994-5 


EWHei 


19 


Kuschnig etal. (1999) 


1995 


RVH^ 


19 


Kuschnig et al. (1999) 


1995 


EWSin 


19 


Kuschnig et al. (1999) 


1995 


EWH(5 


18 


Kuschnig et al. (1999) 


1995 


BeB 


21 


Pyper etal. (1998) 


1998-09 


radio 


5 


Ravi etal. (2010) 


2002-9 


V 


353 


Pojmanski et al. (2001) 


2010-11 


BV 


374 


Henry - this paper 


2010-11 


EWSin 


59 


lanik - this paper 


2010 


V 


38 


Liska - this paper 


2011 


vb 


251 


Janik - this paper 


2011 


UBV 


402 


Janlk & Liska - this paper 



The data used for the analysis of CU Vir are given in Table B.l. 
Here we used the following abbreviations: EW - the equivalent 
width, RV - radial velocity, /3 - Hfi photometry, UBV - Johnson 
UBV photometry, uvby - Stroemgren uvby photometry, ByVtHp 
- Hipparcos photometry, B^tf - the mean longitudinal magnetic 
induction, spf - magnitudes derived from spectrograms obtained 
by UV satellites 0A02 and lUE, and radio - timings of ra- 
diopulses. 



